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I INTRODUCTION 

The Green functions of partial differential operators are of great importance in mathe- 
matical physics and quantum theory |jl|, ^, ^. In particular, the singularities of the 
Green functions on the diagonal play a crucial role in the renormalization procedure of 
quantum field theory 0. 

In this paper we study a special class of differential operators of higher orders that 
appear, in particular, in higher-derivative field theories (for example, in higher-derivative 
quantum gravity 0). Namely, we consider differential operators of higher even orders, 
2N, of the form 

H = {FN + m^)---{F2 + m^){F,+m^), (1) 

where Fi, (z = 1, . . . , N) are Laplace type operators, i.e. second-order operators with the 
scalar leading symbol determined by the metric of the manifold 

= leP, (2) 

where = g'^'^{x)^fi^u, and m is some constant. In the next section we give precise 

definition of Laplace type operators and describe their properties. Obviously, the operator 
H also has a scalar leading symbol 

^l{H) = (3) 

and is a particular case of the so called minimal higher-derivative operators. 

We would like to stress from the beginning that the operators Fi, even if of the same 
Laplace type, are different and, in general, do not commute. They only have the same 
leading, second-derivatives, part, but differ in the first and zero order parts. 

The Green function, or the resolvent kernel, Gh-\{x,x') of the operator H is the 
kernel of the Green, or resolvent, operator Gh-\ = {H — A)^^ and is defined by requring 
it to satisfy the differential equation 

{H-\)Gh-x{x,x') = 6{x,x'), (4) 

with H acting on the first argument of the Green function and x') being the covariant 
delta-distribution. 

In the present paper we will consider mainly the case of compact complete manifolds 
without boundary and assume the metric of the manifold to be positive definite. Then 
the Laplace type operators Fi as well as the operator H are elliptic and eq. @ defines a 
unique Green function. We will also assume for simplicity the constant m to be sufficiently 
large so that all operators {Fi + m^) and also H are positive. 

However, we would like to make some remarks concerning this subject: 

i) First, for manifolds with boundary (in elliptic case) one has to impose additionally 
some suitable boundary conditions to make the solution of the eq. (^ unique. How- 
ever, the singularities of the Green function near the diagonal and the asymptotic 
expansions of it as m ^ oo do not depend on the boundary conditions. 
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ii) Second, when the metric of the manifold has Minkowski type signature, sign (7 = 
(—,+,...,+), the Laplace type operators Fj are hyperbolic. In this case the Green 
function can be fixed by the Wick rotation, i.e. the analytic continuation from the 
Euclidean case, which is equivalent to adding an infinitesimal negative imaginary 
part to the constant m^, — ie. In other words, we consider the Feynman 

propagators. The formulas for the hyperbolic case can be obtained just by analytic 
continuation. 

II LAPLACE TYPE OPERATORS 

Let (M, g) be a smooth Riemannian manifold of dimension d with a metric g. To simplify 
the exposition we assume the manifold M to be complete and compact, i.e. without 
boundary, dM = 0, and the metric g to be positive definite. 

Let V{M) be a smooth vector bundle over the manifold M, End (V) be the bundle of 
all smooth endomorphisms of the vector bundle V, and C°°{M, V) and C°°(M, End (V)) 
be the spaces of all smooth sections of the vector bundles V and End (V^). We assume, as 
usual, the vector bundle V to be Hermitian, i.e. there is a Hermitian pointwise fibre scalar 
product. Then the dual vector bundle V* is naturally identified with V and a natural 

inner product is defined using the invariant Riemannian volume element dvol (x) on 
the manifold M. The completion of C°°{M, V) in this norm defines the Hilbert space of 
square integrable sections of the vector bundle L'^{M, V). 

Let, further, V be a connection, or covariant derivative, on the vector bundle V which 
is compatible with the Hermitian metric on the vector bundle V. Denoting by T*M the 
cotangent bundle we define the tensor product connection on the tensor product bundle 
T*M (g) y by means of the Levi-Civita connection. Similarly, we extend the connection 
V with the help of the Levi-Civita connection to C°°{M,V)-yahied tensors of all orders 
and denote it just by V. 

Usually there is no ambiguity and the precise meaning of the covariant derivative is 
always clear from the nature of the object it is acting on. 

Let, further, tr^ denote the contraction of sections of the bundle T*M^T*M^V with 
the metric on the cotangent bundle, and Q G C°°{M, End{V)) be a a smooth Hermitian 
section of the endomorphism bundle. Then we define the generalized Laplacian 

□ = trgVV (5) 

and a Laplace type differential operator F : C°°{M, V) — > C°°{M, V) by 

F^-n+Q. (6) 

Let x^, {/J, = 1,2, ... ,d), he a. system of local coordinates and and dx^ be the local 
coordinate frames for the tangent and the cotangent bundles. We adopt the notation that 
the Greek indices label the tensor components with respect to local coordinate frame and 
range from 1 through d. Besides, a summation is always carried out over repeated indices. 
Let gi^i, — [d^, d^) be the metric on the tangent bundle, g^'^ — {dx'^, dx^) be the metric 
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on the cotangent bundle, g = det gfj_^, and Afj, be the connection 1— form on the vector 
bundle V. Then the generalized Laplacian reads 

□ = g^'-V.V. = g-"\d, + A,)g"'g^''{d. + A). (7) 

It is worth noting that every second-order differential operator with a scalar 

leading symbol given by the metric tensor is of Laplace type and can be put in the 
form d^) by choosing the appropriate connection 1-form A and the endomorphism Q. 

The operator F has a positive leading symbol and is elliptic. It is symmetric with 
respect to the natural inner product. Moreover, the operator F is essentially self- 
adjoint, i.e. there is a unique 

self-adjoint extension F of the operator F. However, we will not be very careful about 
distinguishing between the operator F and its closure F, and will simply say that the 
operator F is elliptic and self-adjoint. 

II. 1 Heat kernel of a Laplace type operator 

Hence, the operator Up^t) = exp{—tF) for t > is well defined as a bounded operator 
on the Hilbert space L'^{M,V) of square integrable sections of the vector bundle V. 
This operator forms a one-parameter semi-group. The kernel UF{t\x,x') of this operator 
satisfies the heat equation and is called the heat kernel. For t > the heat kernel is a 
smooth section of the external tensor product of the vector bundles V and V* over the 
tensor product manifold M x M. 

In the following we are going to study the Green function and the heat kernel only 
locally (like in |^), i.e. in the neighbourhood of the diagonal of M x M. We will keep a 
point x' of the manifold fixed and consider a small geodesic ball, i.e. a small neighbourhood 
of the point x': B^i = {x G M\r{x,x') < e}, r{x,x') being the geodesic distance between 
the points x and x' . We will take the radius of the ball sufficiently small, so that each 
point X of the ball of this neighbourhood can be connected by a unique geodesic with the 
point x' . This can be always done if the size of the ball is smaller than the injectivity 
radius of the manifold at a;', e < rinj(x'). 

Let cr(x, x') be the geodetic interval, also called world function, defined as one half the 
square of the length of the geodesic connecting the points x and x' 

a{x,x') = -r^{x,x'). (8) 

The first derivatives of this function with respect to x and x' define the tangent vector 
fields to the geodesic at the end points x and x' and the determinant of the mixed second 
derivatives defines a so called Van Vleck-Morette determinant 0, ||, || 

A(a;, x') = g-'/\x)det (-V^V;,a(x, x'))g-'/\x'). (9) 

Let, finally, V{x,x') denote the paralell transport operator along the geodesic from the 
point x' to the point x. It is a section of the external tensor product of the vector bundle 
V and V* over M x M, or, in other words, it is an endomorphism from the fiber of V over 
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x' to the fiber of V over x. Near the diagonal of M x M all these two-point functions are 
smooth single-valued functions of the coordinates of the points x and x' . To simplify the 
consideration, one can assume that these functions are analytic in the ball Bx'- 

The heat kernel of a Laplace type operator is well described by factorizing out the 
semi-classical factor 0, ||, 0] 

Vp{t\x, x') = {A7it)-'^/^A^/\x, x') exp (-^^(a^, a;')) ^F{t\x, x'). (10) 

The function QF{t\x,x') is called the transfer function of the operator F. Obviously, 

[/^+^2(t) = e-'^^'Upit), QF+m-it) = e-'^'npit). (11) 

It can be proved that for a positive operator F + the transfer function Qp^^n'^lt) 
can be presented in form of an inverse Mellin transform (we slightly change notation here 
in comparison to |]^) 

^ c+ioo 

^F+m^{t) = — I dqt'^T{-q)bF+m4q), (12) 



C—lOO 



where c is a negative constant. The function bp^,^2[q) is an entire function of q, so that 
the integrand has simple poles at the positive integer points q = k, [k = 0,1, . . .). The 
function bp+jn'2{q\x,x') satisfies a functional-differential equation 

(^^ + ^DjbF+m4q) = Lbp+m^iq-l), (13) 

where 

L = A-'/^FA'/^ + = ^ ^1/2 + g + ^^4) 

D = a^'V/,, = V^a, (15) 

and the initial condition 

bF+m^{0\x,x') =V{x,x). (16) 

This equation (together with the condition of analyticity and some asymptotic condition 
at g — s> const ± ioo) enables one to compute bp^rn'^{q) if one fixes its initial value at 
some arbitrary point q = go- (For more details, see 0). The initial condition at the 
origin produces the values of bp^^'^lk) at positive integer points, the initial condition at 
the point g = 1/2 would give the values of bp^rn?[q) at all half- integer positive points 
bpj^rn?{k + 1/2), etc. Differentiating the eq. (|13D with respect to g one obtains another 
recursion 



+ ^F+m^io) =Lb'p^^2{q - 1) + ^Dbp+rn^{q), 



(17) 



where 

d 

b'p+m^iq) = ^bp+m^{q). 
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which enables one to compute the derivatives of the function hF+m?{q) at positive integer 
points if one fixes its value &^_,_^2(0). 

Moving the contour of integration to the right and taking into account that the residues 
of the gamma-function T{—q) at the points q = k equal {—l)^/k\ one obtains 0] 



k=0 

where 



nF+m^{t) = J2^^ak{F + m^) + — J dqt'^T{-q)bF+m^{q), (19) 



ak{F + m') = hp+mAk) =\-^\ ^F+^^if) ^ . (20) 



dt 



t=o' 



n is an arbitrary positive integer, n > 1, and c„ ranges in the interval n — 1 < Cn < n. 
The coefficients ak{F + m^) are the famous Hadamard-Minakshisundaram-De Witt-Seeley 
(HMDS)-coefficients [|l], 0, g |Tl[ to the operator F + m^. They can be computed in form 
of covariant Taylor series from the recursion relations that are obtained from the equation 
(p!3|) (by putting q ^ k) with the initial condition ao(-F + m^) = V. The diagonal values 
of the HMDS-coefficients are known in general case up to 04 ||12|, |13|, |14|. In some 
particular cases, e.g. in flat space, there are results for higher-order coefficients. For a 
review of the methods for calculation of HMDS-coefficients 

and further references see |T^ . 

Taking into account eq. (|ll]) we find from (pO]) 

ak{F + m') = X: f^')m'^''"^«"(^)- (21) 



n=0 



Moreover, one can obtain an asymptotic expansion of the function bp+m^iQ) as m — > 00 



W(.)-E „,r^^J-;|,) .»-'-"V,(F). (22) 

From (|T9[) there follows an asymptotic expansion of the transfer function as t ^ 
(called Schwinger-De Witt expansion in the physical literature) [@, |[ |, ^ H 

^F+m<t) ~ + = e"""' E (23) 

k>0 '^^ k>0 

Using this equation one obtains from (0) the well known formula for the asymptotic 
expansion of the functional trace of the heat kernel 

Tr ^2 exp(-tF) ~ {A^)-'IH-'I' ^-^MF), (24) 

A;>0 

where 

AkiF) = Tr L^akiF) = J rfvol (x) tr vak{F\x, x), (25) 



M 



with tr V being the fiber trace. 
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II. 2 Green function of a Laplace type operator 

The Laplace type operator F is a self-adjoint elliptic operator with a real spectrum 
bounded from below. Thus for sufficiently large m the operator {F + m?) is positive 
and its Green operator, or the resolvent operator, Gpj^rn? = {F + rn?)~^ is a bounded 
operator. It can be defined by 



oo 



= J dte-'"" exp(-tF). (26) 



The Green function of the operator {F + m^) is obtained by the kernel form of 
this equation 

oo 

Gf+^2 = (47r)-'^/2^i/2 J dtt-'^/^exp (^~tm^ - ^) fi^ (t). (27) 



This integral converges at the infinity for sufficiently large m. It also converges at 
t = outside the diagonal, i.e. for a 7^ 0. By using the ansatz (|IUp and ([T2|) for the heat 
kernel one obtains a corresponding ansatz for the Green function 



Gp+m^ = (An)-''/' A'/'— J dq (^-j r(-g)r(-g - 1 + d/2)bp+m^{q), 



:28) 



where c < —1/2. 

This ansatz is especially useful for studying the singularities of the Green function, or 
more general, for constructing the Green function as a power series in a. The integrand 
in (|28D is again a meromorphic function. However, contrary to (|12D, we have now a more 
complicated structure of the poles. There are always poles at the points q = k and 
q = k — 1 + d/2, {k = 0, 1, 2, . . .). Here one has to distinguish between odd and even 
dimensions. In odd dimensions, the poles are at the points q = k and q = k + [d/2] — 1/2 
and are simple, whereas in even dimension there are simple poles at g = 0, 1, 2, . . . , (i/2 — 2 
and double poles at the points q = k + d/2 — 1. 

Moving the contour of integration in (^) to the right one can obtain an expansion of 
the Green function in powers of a (Hadamard series). Generally, we obtain 

/^sing I ^non— anal i /^reg fnn\ 

Here G^^^^2 is the singular part which is polynomial in the inverse powers of -/a 

G%'m^ = (47r)-^/^AV2 ^:Ar(d/2 - ^ - 1) ( - ) a,{F + m% (30) 

For the rest we get in odd dimensions 

y^non— anal _i_ /^reg 
^F+m^ -r ^F+m? 
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^1 "-('^+l)/2 „ /^x fc+1/2 

^ ^ ^ ^ fro r(fc + ^)r(A: + |) V2; ^^^^ ^ 



n-(d+l)/2 , , i. 



1 C //t\ <?+l^°'/2 

+ (47r)-^/W/2— y dq{^-^^ r(-g)r(-g-l + rf/2)6^+„.(g), 

(31) 

where ra— l<c„<n — 1/2 and n > id — l)/2. Thus, by putting n ^ oo we recover 
herefrom the Hadamard power series in o for odd d = 1,3,5,... 



~ (-l)'-')'^(4,)-'/^A'/= t - 1 + m. (33) 

In eiJen dimensions, the point is more subtle due to the presence of double poles. 
Moving the contour in ( ]28|) to the right and calculating the contribution of the residues 
at the simple and double poles we obtain 

/^non— anal _i_ /^reg 

H-l)-/-(4.)-/'^A"^.og(^^^)|^^^(|)'a._«,.(f + ™^) 

+ (_lW2-l(4^)-rf/2^1/2y \ ('^y 

^ ^ ^ ^ tok\^{k + d/2)\2) 

X {6i.+^2(A; - 1 + d/2) - [log/x^ + ^(A; + 1) + ^(A; + d/2)] a^_^^,{F + m^)} 
+ (4^)-./2^i/2_ j r(-g)r(-g-l + rf/2)6^+„.(g), 

(34) 

where /i is an arbitrary mass parameter introduced to preserve dimensions, n — 1 < c„ < n 
and '^{z) = {d/dz)log T{z). If we let n ^ cxo we obtain the Hadamard expansion of the 
Green function for even d = 2,4, . . . 

^ (-l)^'-M4.)-'^AV^.og {!^) f (I)' + 

(35) 
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log/i' + ^{k + 1) + ^{k + d/2)] a^_^^{F + m2)| 

(36) 

Note that the singular part (which is a polynomial in inverse powers of and the 
non- analytical parts (proportional to ^/a and log a) are expressed uniquely only in terms 
of the local HMDS- coefficients ak{F + m?), whereas the regular analytical part contains 
the values of the function 6^+^2(9) at half-integer positive points and the derivatives of 
the function hF+m?{.l) at integer positive points, which are not expressible in terms of the 
local information. These objects are global and cannot be expressed further in terms of 
the local HMDS-coefficients. However, they can be computed from the eqs. (|T^) and 
(p!7|) almost by the same algorithm as the HMDS-coefficients in terms of the value of the 
function 6(g) at some fixed go (see 0). 

The regular part of the Green function has a well defined diagonal value and the 
functional trace. It reads in odd dimensions {d = 1, 3, 5, . . .): 

Tr L2 G-l^. = {-lY'^'^^'W-'/'^^^BF^^.{d/2 - 1) 

d~2—2k 

(-l).-../^(4,)-/Vg^i^^^.(F) (37) 
and in even dimensions {d = 2, 4, 6, . . .) 

log + m{d/2) - (T] Arf/2_i(F + m2)| 



where 



(-l)^/-'(4.)-/=| E ^fi^lr- *(./2 - + log '-^^MF) 

+ E m'-^-^'Vik + 1 - d/2)A,{F) , (38) 

k>d/2 ) 

Bp+m^iq) = TiL^bF+m^iq) = J dvo\{x)tTvbF+rn2{x,x), (39) 

M 

and (r= — ^l'(l) = 0.577. . . is the Euler's constant. 
Thus, we see that 

i) all the singularities of the Green function and the non-analytical parts thereof (pro- 
portional to y/a in odd dimensions and to log a in even dimensions) are determined 
by the HMDS-coefficients ak{F); 
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ii) there are no power singularities, i.e. ^^"^^2 = 0, in lower dimensions d = 1,2; 

iii) there is no logarithmic singularity (more generally, no logarithmic part at all) in 
odd dimensions; 

iv) the regular part depends on the values of the function hp^rn^ (g) at half-integer points 
and its derivative at integer points and is a global object that cannot be reduced to 
purely local information like the HMDS-coefficients. 

The logarithmic part of the Green function is very important. On the one hand it 
determines, as usual, the renormalization properties of the regular part of the Green 
function, i.e. the derivative fi{d/dfi)G^p^^2- In particular, 

g for odd d 

On the other hand, it is of crucial importance in studying the Huygens principle. 
Namely, the absence of the logarithmic part of the Green function is a necessary and 
sufficient condition for the validity of the Huygens principle for hyperbolic operators 



I7| , |T8|, [T^. The HMDS-coefficients and, therefore, the logarithmic part of the Green 



function are defined for the hyperbolic operators just by analytic continuation from the 
elliptic case. Thus, the condition of the validity of Huygens principle reads 

E uJn!^(^]^.. (^) ' -^-i^^MF + m^) = 0, (41) 

or, by using 



A;!r(A; + d/2) V2 



By expanding this equation in covariant Taylor series using the methods of one can 
obtain an infinite set of local conditions for validity of the Huygens principle. 

nW(d/2) 

where V is the exterior symmetric tensor product, g is the metric tensor on the tangent 
bundle and < n\ak > denotes the diagonal value of the symmetrized covariant derivative 
of n-th order More explicitly, 

[ad/2-iiF + m^)] = 0, (44) 
[V^arf/2-i(F + m2)] =0, (45) 

[V(^V,)ad/2-i(F + m^)] + ^g^u[ad/2{F + w?)] = 0, (46) 



where the square brackets denote the diagonal value of two-point quantities: [f{x,x')] = 
f{x,x). 
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III HIGHER ORDER OPERATORS 

Let us now describe in detail the class of higher order operators we are going to study. 
We consider again a Riemannian manifold {M,g) and a vector bundle V{M) with a 
connection V. Let A(i), (i = 1, . . . , A^) be a set of different smooth sections of the vector 
bundle T*M (g) End(V^). They define a set of different connections 

V(,) = V + A^i) (47) 

on the vector bundle V{M), and therefore, a set of different generahzed Laplacians 

□i = trgV(i)V(^) =^'^"(V^ + ^(i)^)(V, + ^(i),) 

= □ +^^"(Ai)^V, + V.Aii)^) + g'"'Aii)f,Aii)u. (48) 

Let, further, Qi be a set of different smooth sections of the vector bundle End(y). Then, 
we define a set of different Laplace type operators 

Fi^-Di + Qi (49) 

and a higher order operator of a special form 

= (F^ + m')---(F2 + m')(Fi + m'). (50) 



III.l Algebraical reduction to Laplace type operators 

We are going to study the Green function Gh-x{x, x') of the higher order operator H — \ 
where A is an arbitrary sufficiently 

large negative constant. To do this, we show, first, that the Green operator Gh-\ = 
[H — A)^^ can be reduced to the Green operator of an auxiliary Laplace type operator. 
This will mean that the Green function of the higher order differential operator will be 
reduced to the Green function of a second- order differential operator of Laplace type, 
which is described in previous section. 

Theorem 1 Let F^, {i = 1,...,N), be some positive elliptic operators and H he an 
operator defined by 

H ^{FN + rri')---{F2 + rrt'){Fi + m^), (51) 

where m is a sufficiently large constant. Let In-i be the unit {N — 1) x {N — 1) matrix 
and On-1 be the square {N — 1) x [N — 1) null matrix. Let P and H be square N x N 
matrices defined by the following block form 
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where X is a constant, and let F = Fi 
the operators F^ 



Fi ... 



Fn be a N X N diagonal matrix formed by 
\ 

(53) 



... F, 



N 



Let T be an operator defined by 



T = F -P. 



(54) 

Let Gh~^\ and Gjr+m^ be the Green operators of the operators [H — A) and (JF + m^) 
satisfying the equations 

(H - X)Gh-x = 1 



and 

Then there holds 



( J" + m^)G'jF+„2 = 1. 



Gh-x = trHG 



Proof. Let us introduce some auxiliary operators 

Zi = Gh-x 

Z2 = iF, + m^)Zi 

Zn = {FN-i + m^)ZN-i. 
Then the equation ( |55D can be rewritten as 

{FN + m^)ZN-\GH-x = l. 
Let us collect the operators Zi in a column- vector 



^ Z^ 



V / 



(55) 

(56) 

(57) 
(58) 



(59) 



(60) 



(61) 



The Green operator Gh-x is the first component of the vector Z and can be obtained by 
multiplying with the row- vector nl = (1,0, ...,0) 



Gh-x — ^{Z. 



(62) 
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Further, let us define another column-vector 



n 



N 







(63) 



so that 



n = n^®nl. (64) 

It is not difficlut to show that the equations ( |5UD can be rewritten in a matrix form 

(65) 



Using the Green operator Gjr+mZ of the operator (JF + m^) we obtain 

Finally, from ( p^ ) by taking into account ( p^ ) we get the Green operator Gh-\ 

which proves the eq. (|57|). To prove eq. ( |5BD we note that 



(66) 



(67) 



d_ 
dX 



tiG 



— tr n G^+^a • 



Further, from the definition of the Green operator we have, obviosly, 

d 



dm? 



tr n G^+„2 = -tr n G^^_^ 



Comparing (^) and (p9| ) and using the eq. (p7| ) we convince ourselves that eq. 
correct too. 



(6^ 



(69) 



IS 



Let us make some remarks. First, 'tr' denotes here the usual matrix trace and has 
nothing to do with the trace over the bundle indices which are left intact. Second, the 
operator F has, obviously, the form 



where □ = □ i 
Obviously 



F = -n + Q, (70) 
□ TV and Q = Qi Q) ■ ■ ■ (B Qn have the same diagonal structure. 



□ = trgVV 



(71) 



where V = Vi©---©VArisa new connection which also has the same diagonal structure. 
Therefore, the operator JF reads 



^ = -□ + g - P. 



(72) 
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Remember that the matrix P is constant. Thus, both the operator F and JF are second- 
order elhptic Laplace type differential operators. However, the operator is not self- 
adjoint because the matrix P is not symmetric. In spite of this, the operator (JF + m?) 
for sufficiently large m is non-degenerate, i.e. there exists a well defined Green operator 
Gr+i-n?i which justifies made assumptions. 



III. 1.1 Properties of the matrix P 

The matrix P plays the role of the matrix term and is of great importance in further 
considerations. That is why, we state below some important properties of it. 
First of all, we note that the matrix 11 is nilpotent 

= (73) 

and, of course, degenerate. The matrix P is not degenerate. The matrix P as well as its 
powers belong to the class of two-diagonal matrices. We compute first the powers of this 
matrix. 





' \ 










[X 1 






,---,P"~'= , (74) 



In-2 
Xh 

= XIn. (75) 

It is also not difficult to show that 

det {z-P) = z^ ~ X. (76) 

Therefore, the eigenvalues pj of the matrix P are just the A^-th roots of A 

p, = {-XYl^e^<'^-'y\ j = l,...,iV, (77) 

where |arg(— A)| < vr. The roots Pj lie on the circle \z\ = |A|"^/^. 

The higher powers of the matrix P are expressed in terms of the lower order powers 
(0) and (^ 

piN+j ^ yipj^ (7g) 

where j = 0, 1, 2, . . . , — 1 and i = 0, 1, 2, . . .. Using these formulas we find that for any 
function f{z) that is analytic at the points Pj{X) there holds 



-1 



Cp 
N-1 

= (^N-i-A^-J)P' (79) 
j=o 

where 

1 r z^ 

Cp 
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and the contour Cp goes counter-clockwise around the points pj in such a way that it does 
not contain the singularities of the function f{z) — inside the contour Cp there must be 
no other singularities except for the points pj. Evidently, the coefficients /3„(A; /) possess 
the following 'periodicity' property 

f3n+kN{\ f) = \'Pn{\; f), k = 0,1,2,..., (81) 

SO that there are essentially only N independent coefficients Po, Pi, . . . ,Pn-i- 

The eqs. (^9]) and ( pS]) determine any function of the matrix P. In particular, using 
the traces of the powers of the matrix P 

NX' for n = iN, z = 0, 1, 2, . . . , 
tr = <! (82) 

otherwise 
we can compute the trace of any function of P 

tr/(P)=iV/5^_i(A;/). (83) 

More generally, 

trP"/(P) = iV/5;v-i+n(A;/). (84) 
Herefrom we obtain the coefficients /3„(A; /) 

1 ^ 

1 N 

= _A(_A)(«-^+i)/^^e^-(2^-i)(-+i)/^/(p^.) (85) 
^ i=i 

From (|82|) , it follows, in particular, 

N 

J2p1 = ^ forn ^ 0, ±A^, ±2iV, . . . . (86) 
i=i 



Moreover, using the definition (|5^) we calculate the traces of the product of the powers 
of the matrix P with the matrix 11 

1 d ^1 f A* forn = 2A^ + A^-l, i = 0,l,2,..., 
trnP" = -4-^trP"+i= ''' '. (87) 

[n + L)dA 1^ otherwise 

By using the equation 

|^tr/(P)=tr l^r(P) (88) 

where f'{z) = {d/dz)f{z), and the eq. (0) or using directly the eq. ( ^71) we calculate 
the trace of the product of the function /(P) with the matrix 11 

1 ^ 

tr n/(P) =/3o(A;/) = - (-A)-(^"i)/^Ee^'^^''"'^^''/fe)- (§9) 
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More generally, one can obtain also 



tr Ii{z - P)-' = -^log det {z-P) = (90) 



as well as 

1 /• p*^ 1 ^ 

^TT% J Z /A ly ■ -J 



Cp 



In the next section we will need also the following lemma, which is actually a gener- 
alization of the eq. 



Lemma 1 Let P and 11 he the N x N matrices defined by (p^ ) and Bi, . . . ,Bj\i_2 be 
arbitrary diagonal matrices. Then for any k < N — 2 and arbitrary non-negative integers 
ni, . . . , rife satisfying the condition ni + ■ ■ ■ + Uk < N — 2 there holds 

tr n EiP^^^aP"^ ■ ■ ■ PfcP"^ = 0. (92) 
III. 1.2 Properties of the symbol cr(J^) 

Let us consider now the symbol cr(JF) of the operator JF. It has the form 

a{J^) = a{F) - P, (93) 

where (t{F) is a diagonal matrix. 

That is why let us consider the following matrix 

M = fx-P, (94) 

where fi = fii (B ■ ■ ■ Q) is a. diagonal matrix with some constants fij on the diagonal. 
Note that the matrices fl and P do not commute with each other [jl,P] 7^ 0. 

We are interested in the eigenvalues qj of the matrix M. It is not difficult to calculate 

N 

trM" = tr/i" = 5]/ij", {n = l,2,...,N -1) (95) 

N 

tiM^ = J2fif + {-lfNX, (96) 

N 

det {z-M) = l[{z- fij) - (-l)^A. (97) 

Thus the eigenvalues qj of the matrix M are determined either from one algebraical 
equation of A^-th order 

(/ii - ;z) ■ ■ ■ (/iTV - ^) - A = (98) 
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or from the system of N equations 



TV 



qi-\ I-Qn = (99) 

N 

gf-^ + --- + 4-^ = E/^f-^ (100) 

N 

?f + ■■■ + < = Ef^f + i-^fNX. (101) 



This enables one to calculate the traces 

AT 

tr/(M) = 5:/(g,). (102) 

i=i 

One can get also analogous to ( pOf ) integral representation 

tr f(M) = — I dzf{z)- , 

where Cm is a counter-clockwise contour around the spectrum of the matrix M and 

N N 

A{z,^^,X)=Y. n i^-^k) (104) 

j = l fc = l; k^j 

is a polynomial of order (A^ — 1)- 
Further, by using the equation 

|-tr/(M)=tr ^/'(M), (105) 

where f\z) = {d/dz)f{z), and the obvious equation 

^M=-n, (106) 

which follows from the definition of the matrix M and the eq. ( p2|) , we calculate the trace 
of the product of the function f'{M) with the matrix 11 

trn/'(M) = -^tr/(M). (107) 

Using the eqs. ( pSj ) and (^) we get, in particular 

1 d 

trHM" = — trM"+^ = forn = 0, 1, . . . , - 2 (108) 

n + loX 

trHM^-i = -1 Atr = (-1)^-^ (109) 
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Note that (|1(J8[) is a particular case of the eq. (021). Further, from (PDI) - (|1(J1|) one obtains 



forn = 0,l,2,...,Ar-2. 



(-1)^. 



Further, we get also 



tr n (z - M)-i = —log det {z - M) = -^-^ . , 

(2;-/ii)---(2-;UAr) - (-l)^^A 



tr HM" 



:-i 



fc-i 



1 



(A; - 1)! \dz J {z-fi,)...{z~ M - (-l)^A 

k-N+l 



2=0 



tr HM^ 



1 



1 



T{k~N + 2) V^y/ ivf^i - 1) ■ • • (y/i7v - 1) - 2/^A 

Note that for k < N — 2 the last equation gives zero in accordance with (|109|) . 
This leads to a very simple and convenient representation 



trn/(M) = (-1)^-1^ J dz 

Cm 



(^-/ii)---(^-/i^)-(-l)^A 



We will need also, in particular, the following formula 



tr He 



-tM 



\N-1 



-tz 



Cm 



- ^i) ■ • ■ (;2 - /iiv) - (-l^A 



1 d 



N 



-tM 



9qj 



From ( |108|) it follows that this function is of order 0{t^ ^) as t 0. Using eq. 
obtain 

^ = -(e n (/^^-^.) 



In particular case A = we have 



A=0 



hence, 



and 



tr He" 



d\ 

-tM 



1 



A=0 



N 



N N 
A=0 j=l k=l;k^j 



(110) 
(111) 

(112) 
(113) 
(114) 



(115) 



116) 
we 



These formulas will be useful when studying the heat kernel for the operator 



(117) 

(118) 
(119) 
(120) 
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III. 2 Green function of the operator H — X 

This shows that the problem is reduced to a Laplace type operator with an additional 
matrix structure. Using the usual heat kernel representation for the Green function Q of 
the Laplace type operator JF we have from eq. (^) 

Corollary 1 

oo 

Gh-x = J rfte-*™'trn exp(-tj^), (121) 



or in the kernel form 

oo 

Gh-x = {ATry'/'A'/' J dtr"/' exp i^-tm' - ^) ^t), (122) 



where 

cj^(t) = tr nfi^(t), (123) 
^lj^(t) = Qjr[t\x,x') is the transfer function of the operator ([5^j. 

Thus we see that there is analogous heat kernel representation for the Green function 
of higher-order operator H — X with a new 'transfer function' ujjr(t). 

III. 2.1 Power of a Laplace type operator 

The simplest case is, of course, the case of equal potential terms Qj and the connections 
V«, i.e. 

= Q, V« = V. (124) 
In this case all operators Fj are equal to each other 

Fj=F = -n +Q, (125) 

and the operator H is the A^-th power of the Laplace type operator F + 

H={F + mY- (126) 

Further, in this simple case we have 

T = F -P (127) 
and the operator F commutes with the matrix P 

[F, P] = 0. (128) 

Therefore, 

exp(-tj^) = exp(tP) exp(-tF). (129) 
and, using eqs. (0) and ( p6D we get 

LuAt) = fimpit), fit) = tr nexp(tP), (130) 
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which is determined by ( PT| ) and (^) 

27iiJ z^-X ' ^ ^ ' 

Cp J=i 

On the other hand, expanding the integrand in powers of A and using the formula (|163|) 
we obtain f(t) in form of a power series 

\kj.Nk 

Thus, the function ujjr(t) is of order 0(t^^^) as t ^ 0. 

In case A = we obtain herefrom the usual heat kernel representation for the inverse 
power of the operator [F + m^) 

1 °° 

Gh = (F + m')-^ = J dtt^-^e"'""' exp{-tF). (133) 







III. 2. 2 Constant perturbations of the potential terms 

Let us consider now the case when the Laplace type operators Fj differ from each other 
only by a constant part of the potential terms, i.e. the connection V'--^-* is the same, but 
the endomorphisms Qj are different 

Qj=Q + fij, (134) 

so that 

= F + /i,-, (135) 

where 

F = -n+Q. (136) 

The operator H has the form 

H={F + m^ + fXN)---{F + m^ + /i2)(F + + /ii). (137) 

The operators Fj, even if different, commute. The corresponding operator has the form 

7 = F + M, (138) 

where M = jl — P and /i = /ii©---©//ivisa diagonal matrix with the constants Hj on 
the diagonal. The matrix M commutes with the operator F 

[F, M] = 0. (139) 

Therefore, 

uj^it) = h{t)nFit), hit) = tr nexp(-tM). (140) 
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Using the eigenvalues of the matrix M, which are determined by the eq. 
this function in the form 



we obtained 



h{t) 



-1 



2m 



-tz 



dz ■ 



Cm 



z-/ii)---(z-M-(-l)^A 



TV 



:i4i) 



As it has been shown in section 3.2, this function is of order 0{t^ ^) as t — > 0. For A = 
we obtained 



h{t) 



N 

= E 

A=0 j=l 



n ^ 



(142) 



III. 2. 3 General case 



As we have seen, in particular cases studied before the function ujjr(t) defined by ( |123| ) is 
of order 0{t^~^) a.s t —>■ 0. It is almost clear that this is valid also in general, since the 
main contribution to the asymptotics give the constant background fields. This can be 
formulated in form of a lemma. 



Lemma 2 Let T he the Laplace type operator defined by ( \54 ) and H be the matrix defined 
by (|5^. Then for its HMDS- coefficients ak{T) there holds 



trnafc(J^)=0 for A; = 0,1,..., A^- 2. 



(143) 



In other words, for the asymptotic expansion as t ^ of its transfer function ujjr[t) there 
holds 

^ E^^c^(•^). (144) 



where 



{N-iy.k\ 



k\ 



{N + k-l)\ 



tr liak+N-ii,^)- 



(145) 



Proof. The coefficients at can be calculated in form of a covariant Taylor series near 
the diagonal (for details and notation see 0) 



n>0 



(146) 



where |n > is the covariant Taylor basis and < n\ak > are the diagonal values of the 
symmetrized covariant derivatives. < n\ak > have the following general form 

<n\ak>= c{ni, . . . ,nk-i) < n\M\nk-i >< nk-i\M\nk-2 > ■ ■ ■ < ni\M\0 >, 

ni,...,nfc_i>0 

(147) 
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where c{ni, . . . ^rik^i) are some constants and < m\M\n > are some End (y)-valued ten- 
sors that vanish for n > m + 2 and for n = m + 1. Important for us is that, since the 
matrix P is constant, it appears only in < n\M\n > 

<n\M\n>= P + B, (148) 

where S is a diagonal matrix. All others tensors < n\M\m > do not depend on the 
matrix P and are diagonal — they are actually polynomials in the curvatures and their 
covariant derivatives 0. So, all the terms in the sum (|147|) have the form 

fiiP"! . . . fi,P"^ j = l,2,...,fc, (149) 

where Bi are some diagonal matrices and the total power of the matrix P is restricted 
hj rii + ■ ■ ■ + rij < k. Therefore, from (^) we find that trll < n\ak >= tr Ilafc = for 
k < N — 2, which proves the lemma. 



Actually, using (|148|) and (^7\) one can prove more, namely, 

co(j^) = tr n ttN-iiJ^) = tr n P^~^ = 1. (150) 

III. 2. 4 Hadamard expansion of the Green function of the operator H — \ 

The behavior of the Green function of Laplace type operators near diagonal is studied in 
the previous section. So, by using the eqs. of sect. 2 we obtain from (^) the Hadamard 
series for the Green function of the higher order operator H — X 

_ rising . ^non-anal , ^reg /-■ r-, \ 

where the singular part reads 



(152) 

Further, we get in odd dimensions: 

^ ' ^ ' {N-iy.to r(A;-iv + i + ^)r(fc + |) ^ ' 



and in even dimensions: 



^^on anal ^ (4^)-<^/2^^1 E ' (155) 
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X {v'^+m-^{k -N + d/2) - [log/i^ + ^(A; + 1) + ^(A; + d/2)] C;^_^+|(J^ + m^)} 

(156) 

where 

CA^ - DTfo + 1) 

^^W(g) = r(g + iV) n6^+„.(g + iV - 1). (157) 

Thus, we find that 

i) the structure of the diagonal singularities of the Green function Gh~\ is determined 
by the coeffcieints Cfc(JF + m^); 

ii) in odd dimenions there is no logarithmic singularity; 

iii) for N > d/2 there are no singularities at all and the Green function is regular on the 
diagonal, i.e. there exists finite coincidence limit G'^'^\ = Gh-x{x, x) = GJ^^ ;^(x, x); 

iv) for = [{d + l)/2] there are no singularities for odd dimension d = 2N — 1 and 
there is only logarithmic singularity in even dimension d = 2N. 

v) The condition for the validity of the Huygens principle for the operator if — A in 
even dimensions reads 

°° 1 /cr\^ 



k 



Similarly to (^31) this also produces an infinite set of local conditions. 



III. 3 Heat kernel of the operator H 

The heat semigroup for the operator H can be obtained by using the standard formula 



exp{-tH) = --^ / dXe-'^Gn^ 

2711 J 



Ch 

_^ C+ioO / \ fc 



L_ f rfAe-*M|-) Gh-x, (159) 



2ni J \d\ 



C—lOO 



where the contour of integration Gh goes counter-clockwise around the spectrum of the 
operator H, k is some sufficiently large integer and c is a negative constant. Since the 
operator ii is a positive self-adjoint operator its spectrum lies on the positive real half- 
axis, i.e. the contour Gh should go from + 00 to —ie + 00 enclosing the origin. Using 
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the Theorem 1 one can express the heat semigroup for the higher order operator H in 
terms of the heat semigroup of the Laplace type operator JF. We have first 

exp{-tH) = — I d\e-'^i — \ trnG^+^2, (160) 

c—ioo 

Now, using the heat kernel representation of the Green operator Gjrj^^2 we obtain 

Corollary 2 The heat semigroup for the operator H is determined by the heat semigroup 
for the operator T 

,-k ^+;°° y / 

exp(-tH) = / dXe^'^ / dse-'""' — tr H exp(-s^) (161) 

27ri J J \oX I 

c-ioo ^ ^ 

and, consequently, the heat kernel for the operator H is given by 

Unit) = (4vr)-'/W/2i- J dXe'^' J ds s-'^' exp [-sm' - uj^s). 



(162) 

We studied the expansion of the Green function of a Laplace type operator in previous 
section in detail. The important fact is, that the HMDS-coefhcients afc(^ + m^) for the 
Laplace type operator (jF + m^) are polynomial in A, they are, of course, also polynomial 
in m^. This means that the singular part of the Green function (including the logarithmic 
singularity), which are expressed in terms of ak{J^ + m^), are also polynomial in A. On 
the other hand, it is well known that 

Ch ^ ^ 

For integer 2; = A; = 0, 1, 2, ... we have, therefore, 

1 

Ch 



d\e'^{-\f = (164) 



Thus, the singular part, which is polynomial in A, does not contribute to the integral 
( 11601) . Therefore, one can substitute in this equation instead of the whole Green function 
only the regular part of it G"''^^, go to the diagonal and take the trace 

TTL.exp{-tH) = — I dAe-*M— j Tr^^trHG^^^,, (165) 

c—ioo 

Finally, using the eqs. (|37D and ([38| ) we have: in odd dimensions 
IV,.expH//)^ - / dAe-y *.wW2-iV). 

C—ioo 

(166) 
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and in even dimensions 



c+ioo 



(_l)d/2-lU ^^k 



C—lOO 




where $^+^2 (g) = Tr (q) and $^+„2(g) = (rf/rfg)$^+^2(g). 

These formulas may be used, for example, to compute the asymptotic expansion of 
the functional trace of the heat kernel for the higher order operator H, which is known 
to have the following general form 

Tr,.exp(-tff) ~ (4vr)-'^/^ ^^yg^^|] r'^^^^^ l^f""^'^^^- ^'^^^ 

We normalize the coefficients Ek{H) so that Eq{H) = Tr^2 1. For = 1 the coefficients 
Ek are determined by the HMDS-coefficients Ek = [{—l)'^/k\]Ak. However, one should 
note that to compute the coefficients Ek{H) one needs to know the asymptotics of the 
function ujjr{t) (or the functions $;F+m2(c?/2 — 1) and $^_,_^2((i/2 — 1)) as A — oo, which 
is unknown, in general, in terms of the known asymptotic expansion as t — 0, i.e. in 
terms of the HMDS-coefficients for a Laplace type operator. The known asymptotics of 
these functions as s — > and m ^ oo do not contribute at all in the eqs. ( |162| ), (|166|) 
and (IT67D. 
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